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ABSTRACT 

Sphere packings in high dimensions interest mathematicians and physicists and have 
direct applications in communications theory. Remarkably, no one has been able to provide 
exponential improvement on a 100-year-old lower bound on the maximal packing density due 
to Minkowski in ci-dimensional Euclidean space W 1 . The asymptotic behavior of this bound 
is controlled by 2~ d in high dimensions. Using an optimization procedure that we introduced 
earlier [TS02I and a conjecture concerning the existence of disordered sphere packings in W d , 
we obtain a conjectural lower bound on the density whose asymptotic behavior is controlled 
by 2~°- 77865 -" d , thus providing the putative exponential improvement of Minkowski's bound. 
The conjecture states that a hard-core nonnegative tempered distribution is a pair correlation 
function of a translationally invariant disordered sphere packing in M. d for asymptotically 
large d if and only if the Fourier transform of the autocovariance function is nonnegative. 
The conjecture is supported by two explicit analytically characterized disordered packings, 
numerical packing constructions in low dimensions, known necessary conditions that only have 
relevance in very low dimensions, and the fact that we can recover the forms of known rigorous 
lower bounds. A byproduct of our approach is an asymptotic conjectural lower bound on the 
average kissing number whose behavior is controlled by 2 - 22134 - , which is to be compared 
to the best known asymptotic lower bound on the individual kissing number of 2 0,2075 - d . 
Interestingly, our optimization procedure is precisely the dual of a primal linear program 
devised by Cohn and Elkies |CE03| to obtain upper bounds on the density, and hence has 
implications for linear programming bounds. This connection proves that our density estimate 
can never exceed the Cohn-Elkies upper bound, regardless of the validity of our conjecture. 



1 Introduction 



A collection of congruent spheres in d- dimensional Euclidean space M. d is called a sphere 
packing P if no two of the spheres have an interior point in common. The packing density or 
simply density 4>(P) of a sphere packing is the fraction of space M. d covered by the spheres. 
We will call 

</>max = SUp (f>(P) (1) 



the maximal density, where the supremum is taken over all packings in M. d . The sphere 
packing problem seeks to answer the following question: Among all packings of congruent 
spheres, what is the maximal packing density maX ; i-e., largest fraction of M. d covered by the 
spheres, and what are the corresponding arrangements of the spheres |Rog64, CS98 ? The 
sphere packing problem is of great fundamental and practical interest, and arises in a variety 
of contexts, including classical ground states of matter in low dimensions |CL95j . the famous 
Kepler conjecture for d = 3 [HalQHj, error-correcting codes |CS98j and spherical codes (CS98 . 

For arbitrary d, the sphere packing problem is notoriously difficult to solve. In the case 
of packings of congruent <i-dimensional spheres, the exact solution is known for the first three 
space dimensions. For d — 1, the answer is trivial because the spheres tile the space so 
that max = 1. In two dimensions, the optimal solution is the triangular lattice arrangement 
(also called the hexagonal packing) with max = n/^/l2. In three dimensions, the Kepler 
conjecture that the face-centered cubic lattice arrangement provides the densest packing with 
0max = tt/V18 was only recently proved by Hales |Hal98j . For 3 < d < 10, the densest known 
packings of congruent spheres are lattice packings (defined below). However, for sufficiently 
large d, lattice packings are likely not to be the densest. Indeed, this paper suggests that 
disordered sphere arrangements may be the densest packings as d — > oo. 

We review some basic definitions that we will subsequently employ. A packing is saturated 
if there is no space available to add another sphere without overlapping the existing particles. 
The set of lattice packings is a subset of the set of sphere packings in M. d . A lattice A in R d is 
a subgroup consisting of the integer linear combinations of vectors that constitute a basis for 
Mr. A lattice packing Pi is one in which the centers of nonoverlapping spheres are located at 
the points of A. In a lattice packing, the space M. d can be geometrically divided into identical 
regions F called fundamental cells, each of which contains the center of just one sphere. Thus, 
the density of a lattice packing <p L consisting of spheres of unit diameter is given by 

,l Ml/2) 



Vol(F) 
2 



(2) 



where 

7T d/2 

Vl (R) = / - R d (3) 
u ; r(l + d/2) v ; 

is the volume of a d- dimensional sphere of radius R and Vol(F) is the volume of a fundamental 

cell. We will call 

0Lx = sup 0(P L ) (4) 

P L CM d 

the maximal density among all lattice packings in M d . For a general packing of spheres of 
unit diameter for which a density 4>(P) exists, it is useful to introduce the number (or center) 
density p defined by 

4>(P) 

which therefore can be interpreted as the average number of sphere centers per unit volume. 

Three distinct categories of packings have been distinguished, depending on their be- 
havior with respect to nonoverlapping geometric constraints and/or externally imposed vir- 
tual displacements: locally jammed, collectively jammed, and strictly jammed |TS01l rTDS03, 
DTSC04J. Loosely speaking, these jamming categories, listed in order of increasing stringency, 
reflect the degree of mechanical stability of the packing. A packing is locally jammed if each 
particle in the system is locally trapped by its neighbors; i.e., it cannot be translated while 
fixing the positions of all other particles. Each sphere simply has to have at least d + 1 con- 
tacts with neighboring spheres, not all in the same d-dimensional hemisphere. A collectively 
jammed packing is any locally jammed configuration in which no finite subset of particles can 
simultaneously be continuously displaced so that its members move out of contact with one 
another and with the remainder set. A strictly jammed packing is any collectively jammed 
configuration that disallows all globally uniform volume-nonincreasing deformations of the 
system boundary. Importantly, the jamming category is generally dependent on the type of 
boundary conditions imposed (e.g., hard wall or periodic boundary conditions) as well as the 
shape of the boundary. The range of possible densities for a given jamming category decreases 
with increasing stringency of the category. Whereas the lowest- density states of collectively 
and strictly jammed packings in two or three dimensions are currently unknown, one can 
achieve locally jammed packings with vanishing density [B6r64 . This classification of pack- 
ings according to jamming criteria is closely related to the concepts of "rigid" and "stable" 
packings found in the mathematics literature CBB98 . 

In the next section, we summarize some previous upper and lower bounds on the maxi- 
mal density. For large d, the asymptotic behavior of the well-known Minkowski lower bound 
|Min05j on the maximal density is controlled by 2~ d . Thus far, no one has been able to provide 



3 



exponential improvement on this lower bound. Using an optimization procedure and a con- 
jecture concerning the existence of disordered sphere packings in high dimensions, we obtain 
conjectural lower bounds that yield the long-sought asymptotic exponential improvement on 
Minkowski's bound. We believe that consideration of truly disordered packings is the key 
notion that will yield exponential improvement on Minkowski's lower bound. A byproduct of 
our approach is an asymptotic conjectural lower bound on the average kissing number that is 
superior to the best known asymptotic lower bound on the individual kissing number. 



2 Some Previous Upper and Lower Bounds 

The nonconstructive lower bounds of Minkowski [Min05| established the existence of reason- 
ably dense lattice packings. He found that the maximal density max among all lattice packings 
for d > 2 satisfies 

6 L > C M ( 6 ) 

where ((d) = Yl'kLi k~ d is the Riemann zeta function. Note that for large values of d, the 
asymptotic behavior of the Minkowski lower bound is controlled by 2~ d . Since 1905, many 
extensions and generalizations of © have been obtained |DR47| lBa!92|. IUS98j . but none of 



these investigations have been able to improve upon the dominant exponential term 2~ d . It 
is useful to note that the density of a saturated packing of congruent spheres in R d for all d 
satisfies 

The proof is trivial. A saturated packing of congruent spheres of unit diameter and density 
in M d has the property that each point in space lies within a unit distance from the center of 
some sphere. Thus, a covering of the space is achieved if each sphere center is encompassed 
by a sphere of unit radius and the density of this covering is 2 d <ft > 1. Thus, the bound (|7|). 
which is sometimes called the "greedy" lower bound, has same the dominant exponential term 
as (JHJ). In Section 4.1, we show that there exists a construction of a disordered packing of 
congruent spheres that realizes the weaker lower bound of (J7J), i.e., <ft — ^ 

The best currently known lower bound on max was obtained by Ball |Bal92j . He found 

that 

L 2(d-l)((d) f . 



As d — > oo, Table [U gives the dominant asymptotic behavior of several lower bounds on 
for large d. 



max 
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Table 1: Dominant asymptotic behavior of lower bounds on 0^ ax for large d. 



(2)2" d 


Minkowski (1905) 


[\n{V2)d}2- d 


Davenport and Rogers (1947) 


{2d)2~ d 


Ball (1992) 



Nontrivial upper bounds on the maximal density max for any sphere packing in R rf 
have been derived. Blichfeldt [B li29j showed that the maximal packing density for all d 
satisfies </> max < (d/2 + l)2~ d ^. This upper bound was improved by Rogers Rog58] |Rog64| 
by an analysis of the Voronoi cells. For large d, Rogers' upper bound asymptotically becomes 
d2~ d l 2 /e. Kabatiansky and Levenshtein |KL78j found an even stronger bound, which in the 
limit d -> oo yields max < 2-°- 599M ( 1+ °( 1 ». Cohn and Elkies |CK03| obtained and computed 
linear programming upper bounds, which provided improvement over Rogers' upper bound for 
dimensions 4 through 36. They also conjectured that their approach could be used to prove 
sharp bounds in 8 and 24 dimensions. Indeed, Cohn and Kumar fUK04| used these techniques 
to prove that the Leech lattice is the unique densest lattice in 3ft 24 . They also proved that no 
sphere packing in 9ft 24 can exceed the density of the Leech lattice by a factor of more than 
1 + 1.65 x 10~ 30 , and gave a new proof that E 8 is the unique densest lattice in 9ft 8 . Tabled 
provides the dominant asymptotic behavior of several upper bounds on max for large d. 

Table 2: Dominant asymptotic behavior of upper bounds on max for large d. 



(d/2)2-°- 5d 


Blichfeldt (1929) 


(d/e)2-°- 5d 


Rogers (1958) 


2-0.5990<i 


Kabatiansky and Levenshtein (1979) 



3 Realizability of Point Processes 

As will be described in Section our new approach to lower bounds on the density of sphere 
packings in M. d rests on whether certain one- and two-point correlation functions are realizable 
by sphere packings. As we will discuss, a sphere packing can be regarded to be a special 
case of a point process and so a more general question concerns the necessary and sufficient 
conditions for the realizability of point processes in M. d . Before discussing the realizability of 
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point processes, it is useful to recall some basic results from the theory of spatial stochastic 
(or random) processes. Let x = (xi, X2, . . . Xd) denote a vector position in M. d . Consider a 
stochastic process {5^(x; u) : x G M. d ; uj G Q}, where K(x; u) is a real-valued random variable, 
w is a realization generated by the stochastic process and (Q, V) is a probability space 
(i.e., Q is a sample space, T is a a-algebra of measurable subsets of Q, and V is a probability 
measure). For simplicity, we will often suppress the variable uj. 

3.1 Ordinary Stochastic Processes 

We will assume that the stochastic process is translationally invariant (i.e., statistically ho- 
mogeneous in space). Let us further assume that the mean \i = (V(x)) and autocovariance 
function 

x(r) = ([y(x)-Ai][r(x + r)-/i]> (9) 

exist, where angular brackets denote an expectation, i.e., an average over all realizations. The 
fact that the mean \i and autocovariance function %(r) are independent of the variable x is a 
consequence of the translational invariance of the measure. Clearly, 

X (0) = (Y 2 ) - f (10) 

and it follows from Schwarz's inequality that 

| X (r)|<(F 2 )-/x 2 . (11) 

It immediately follows |Loe63j that the autocovariance function x{ T ) m ust be positive semidef- 
inite (nonnegative) in the sense that for any finite number of spatial locations ri, r2, . . . , r m 
in M d and arbitrary real numbers a±, 0,2, ■ ■ ■ , a m , 

m m 

^^a^&fr.-r^O. (12) 

*=i 3=1 

It is clear that ([^(x + r) — F(x)] 2 ) = 2[x(0) — x{ r ))- Thus, if the autocovariance function 
x( T ) is continuous at the point r = 0, the process K(x) on M. d will be mean square continuous, 
i.e., lim r ^ ([^( x + r) — ^( x )] 2 ) = for all x. Stochastic processes that are continuous in the 
mean square sense will be called ordinary. It is simple to show that if x( r ) is continuous at 
r = 0, it is continuous for all r. 

Does every continuous positive semidefinite function /(r) correspond to a translationally 
invariant ordinary stochastic process with a continuous autocovariance x( r )? The answer is 
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yes, and a proof is given in the book by Loeve Loe63j for stochastic processes in time. Here 
we state without proof the generalization to stochastic processes in space. 



Theorem 3.1 A continuous function /(r) on M. d is an autocovariance function of a transla- 
tionally invariant ordinary stochastic process if and only if it is positive semide finite. 

Remark: 

Assuming that /(r) is positive semi definite, one needs to show that there exists a random 
variable Y(x) on M. d such that ([Y(x) — ^t][Y(x+r)— fi\) = /(r). This is done by demonstrating 
that there exists a Gaussian (normal) process for every autocovariance function [Loe63]. A 
crucial property of a Gaussian process is that its full probability distribution is completely 
determined by its mean and autocovariance. 

The nonnegativity condition (JP2j) is difficult to check. It turns out that it is easier to es- 
tablish the existence of an autocovariance function by appealing to its spectral representation. 
We denote the space of absolutely integrable functions on M d by L 1 . The Fourier transform 
of an L 1 function / : M. d — > 3? is defined by 

/(k)= / f(r)e- ik - r dv. (13) 

JK. d 

If / : M d — > K. is a radial function, i.e., / depends only on the modulus r = |r| of the vector 
r, then its Fourier transform given by 

/» = (fer)« r^m ^-f!} dr, (14) 
J o I rcr l 

where k is the modulus of the wave vector k and J u (x) is the Bessel function of order v. The 
Wiener-Khint chine theorem states that a necessary and sufficient condition for the existence 
of a continuous autocovariance function %(r) of a translationally invariant stochastic process 
{Y(x) : x G M. d } is that its Fourier transform is nonnegative everywhere, i.e., x(k) > for all 
k |Yag87[ lTo r02j. The key "necessary" part of the proof of this theorem rests on a well-known 
theorem due to Bochner |Boc36j . which states that any continuous function /(r) is positive 
semidefinite in the sense of (fT2|) if and only if it has a Fourier-Stieltjes representation with a 
nonnegative bounded measure. 
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3.2 Generalized Stochastic Processes 



The types of autocovariance functions that we are interested in must allow for generalized 
functions, such as Dirac delta functions. The Weiner-Khintchine theorem has been extended 
to autocovariances in the class of generalized functions called tempered distributions, i.e., 
continuous linear functionals on the space S of infinitely different iable functions $(x) such 
that $(x) as well as all of its derivatives decay faster than polynomially. Nonnegative tempered 
distributions are nonnegative unbounded measures v on R d such that 



for some n. The interested reader is referred to the books by Gel'fand |GV64j and Ya- 



{Y(<&(x)) : x G M d } is a generalized stochastic process if for each $(x) G S we have a random 
variable K(<&(x)) that is linear and mean square continuous in $. Then the mean is the linear 
functional /z($(x)) = (Y($i(x))) and the autocovariance function is the bilinear functional 
([F($i(x)) -/x($(xi))][F($ 2 (x + r)) -/x($(x 2 ))]), which we still denote by x( r ) fc> r simplicity. 

Theorem 3.2 A necessary and sufficient condition for an autocovariance function %(r) to 
come from a translationally invariant generalized stochastic process {Y($(x)) : x G M. d } is 
that that its Fourier transform x(k) is a nonnegative tempered distribution. 

Remark: 

We will call Theorem 13.21 the generalized Wiener-Khintchine theorem. 
3.3 Stochastic Point Processes 

Loosely speaking, a stochastic point process in M. d is defined as a mapping from a probability 
space to configurations of points x 1 ,x 2 ,x 3 ... in IR d . More precisely, let X denote the set of 
configurations such that each configuration x G X is a subset of M. d that satisfies two regularity 
conditions: (i) there are no multiple points (xj ^ Xj if i ^ j) and (ii) each bounded subset 
of lR d must contain only a finite number of points of x. We denote by N(B) the number of 
points within xdB, B G B, where B is the ring of bounded Borel sets in R d . Thus, we always 
have N(B) < oo for B G B but the possibility A^(M a! ) = oo is not excluded. We denote by 
U the minimal u-algebra of subsets of X that renders all of the functions N(B) measurable. 




(15) 



glom |Yag87j for details about generalized stochastic processes. It suffices to say here that 
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Let (O, J 7 , V) be a probability space. Any measurable map x(oj) : Q — > X, uj G Q, is called a 
stochastic point process |SKM95j . Point processes belong to the class of generalized stochastic 
processes. 

A particular realization of a point process in R d can formally be characterized by the 
random variable 

oo 

n(r) = ^5(r-x J ) (16) 

i=l 

called the "local" density at position r, where S(r) is a ^-dimensional Dirac delta function. 
The "randomness" of the point process enters through the positions x 1; x 2 , . . .. Let us call 

\\% (17) 

the indicator function of a measurable set A C R d , which we also call a "window." For a 
particular realization, the number of points N(A) within such a window is given by 



N(A) = / n(r)I A (r)dr 

oo „ 

= V / 5(r-Xi)I A (r)dr 



i=l 
i>l 

Note that this random setting is quite general. It incorporates cases in which the location of 
the points are deterministically known, such as a lattice. A packing of congruent spheres of 
unit diameter is simply a point process in which any pair of points cannot be closer than a 
unit distance from one another. 

It is known that the probability measure on (X,U) exists provided that the infinite set 
of n-point correlation functions p n , n = 1,2,3... meet certain conditions |Len73| ILen75at 
ILen75bj . The n-point correlation function p n (i"i, . . . , r n ) is the contribution to the ex- 
pectation (n(r 1 )n(r 2 ) ■ ■ -n(r n )) when the indices on the sums are not equal to one another, 
i.e., 

p n (ri,r 2 , ...,r n ) = ( ^ <J(ri - x ix )5(r 2 - x i2 ) • • • S(r n - x in ) \ . (19) 

Note that the distribution-valued function p n (r 1; r 2 , . . . , r n ) also has a probabilistic interpreta- 
tion: apart from trivial constants, it is the probability density function associated with finding 
n different points at positions i"i,r 2 , . . . ,r n . For this reason, p n is also called the n-particle 
density and, for any n, has the nonnegativity property 

p n (n,r 2 ,...,r n ) > Vr,GR d (i = 1, 2, . . .n). (20) 
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Translational invariance means that for every constant vector y in M. d , p n (r 1 ,r 2 , . . . ,r n ) = 
Pn( r i + y ? • • • , r n + y), which implies that 

p n (ri, r 2 , . . . , r n ) = p n g„{r 12 , r ln ), (21) 

where p is the number (or center) density and g„(i"i 2 , . . . , ri„) is the n-particle correlation 
function, which depends on the relative positions ri 2 , ri 3 , . . ., where = rj — and we have 
chosen the origin to be at r±. 

For such point processes without long-range order, g n (r 12 , . . . , r ln ) — > 1 when the points 
(or "particles") are mutually far from one another, i.e., as |r^-| — > oo (1 < % < j < oo), 
p n ( r i) r 2, • • • , r n) p n ■ Thus, the deviation of g n from unity provides a measure of the degree 
of spatial correlation between the particles, with unity corresponding to no spatial correlation. 
Note that for a translationally invariant Poisson point process, g n is unity for all values of its 
argument. 

As we indicated in the beginning of this section, the first two correlation functions, 
/0i(i"i) = p and p 2 ( r i, r 2) = P 2 fi f 2( r ) ) for translationally invariant point processes are of central 
concern in this paper. If the point process is also rotationally invariant (statistically isotropic), 
then g 2 depends on the radial distance r = |r| only, i.e., 

<7 2 (r) = <7 2 (r), (22) 

and is referred to as the radial distribution function. Strictly speaking, one should use different 
notation for the left and right members of Eq. (22), but to conform to conventional statistical- 
mechanical usage, we invoke the common notation for both. Because p 2 (ri,r 2 )/p = pg 2 (r) is 
a conditional joint probability density, then 

Z(r 1 ,r 2 )= psi(r)g 2 (r) dr (23) 

Jri 

is the expected number of points at radial distances between r\ and r 2 from a randomly chosen 
point. Here S\{f) is the surface area of a rf-dimensional sphere of radius r given by 

2 n d/2 r .d-i 

siM = tw^t (24) 

For a packing of congruent spheres of unit diameter, g(r) — for < r < 1, i.e., 

snpp(g 2 ) C {r : r > 1}. (25) 

Note that the radial distribution function g 2 (r) (or any of the p n ) for a point process must be 
able to incorporate Dirac delta functions. We will specifically consider those radial distribution 
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functions that are nonnegative distributions. For example, g 2 {r) f° r a lattice packing is the 
rotational symmetrization of the sum of delta functions at lattice points at a radial distance 
r from any lattice point [TS03 . 

For a translationally invariant point process, the autocovariance function x(r) is related 
to the pair correlation function via 

X (r) = p5(r)+p 2 h(r), (26) 

where 

h(r) = g 2 (r) - 1, (27) 

is the total correlation function. This relation is obtained using definitions (jHJ) and (jlfij) with 
F(x) = n(x). Note that x( r ) — P<K r ) (i-e., h = 0) for a translationally invariant Poisson point 
process. Positive and negative pair correlations are defined as those instances in which h is 
positive (i.e., g% > 1) and h is negative (i.e., g 2 < 1), respectively. The Fourier transform of 
the distribution- valued function x(r) is given by 

x(k)=p + p 2 Mk), (28) 

where h(k) is the Fourier transform of h(r). It is common practice in statistical physics to 
deal with a function trivially related to the spectral function x(k) called the structure factor 
S(k), i.e., 

S(k) = ^ = l + p/i(k). (29) 

A natural question to ask at this point is the following: Given a positive number density p 
and a pair correlation function (72(1*), does there exist a translationally invariant point process 
in M. d with measure V for which p and g 2 are one-point and two-point correlation functions, 
respectively. Two obvious nonnegativity conditions |TS02j that must be satisfied are the 
following: 

g 2 (r) > for all r (30) 

and 

S(k) = 1 + ph(k) > for all k. (31) 

The first condition is trivial and comes from (|2U|) with n = 2. The second condition is 
nontrivial and derives from the generalized Wiener-Khintchine theorem ()3.2j) using relations 
(|28|) and (J29|) . However, for realizability of point processes in arbitrary dimension d, the 
two standard conditions ()30|) and (J3~T|) are only necessary, not necessary and sufficient. The 
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same state of affairs applies to the theory of random sets |Tor02j . where it is known that the 
Wiener- Khint chine theorem only provides a necessary condition on realizable autocovariance 
functions. The simplest example of a random set is one in which Mr is partitioned into two 
disjoint regions (phases) but with an interface that is only known in a probabilistic sense. (A 
packing can therefore be viewed as a special random set). Thus, a random set is described 
by a random variable that is the indicator function for a particular phase, i.e., it is a binary 
stochastic process. The class of autocovariances that comes from a binary stochastic process 
is a subclass of the total class that comes from an ordinary process {Y(x) : x G M d } and 
meets the existence condition of Theorem 13.21 Similarly, the class of autocovariances that 
comes from a point process is a subset of of a generalized process {Y($(x)) : x G M d } and 
therefore the existence condition of Theorem 13.21 is only necessary. 

It has recently come to light that a positive g-i for a positive p must satisfy an uncount- 
able number of necessary and sufficient conditions for it to correspond to a realizable point 
process [CL04 . However, these conditions are very difficult (or, more likely, impossible) to 
check for arbitrary dimension. In other words, given p\ = p and pi = p 2 g2, it is difficult to 
ascertain if there are some higher-order functions p 3 , p 4 , . . . for which these one- and two-point 
correlation functions hold. Thus, an important practical problem becomes the determination 
of a manageable number of necessary conditions that can be readily checked. 

One such additional necessary condition, obtained by Yamada |Yam61j , is concerned with 
<7 2 (A) = ((N(A) — (N(A))) 2 ), the variance in the number of points N(A) contained within a 
window A C M d . Specifically, he showed that 

a 2 (A) > 0(1 -0), (32) 

where is the fractional part of the expected number of points p\A\ contained in the window. 
This inequality is a consequence of the fact that the number of points N(A) within a window 
at some fixed position is an integer, not a continuous variable, and sets a lower limit on the 
number variance. We note in passing that the determination of the number variance for lattice 
point patterns is an outstanding problem in number theory |Ken48t lKR53| l5"S05j . The number 
variance for a specific choice of A is necessarily a positive number and generally related to the 
total pair correlation function h(r) for a translationally invariant point process [TS03J. In the 
special case of a spherical window of radius R in M d , it is explicitly given by 



a 2 {R) = pv 1 {R) 



1+p /i(r)«2(r; R) dr 



> 0(1-0), (33) 



where u 2 (R) is the number variance for a spherical window of radius R in R d , vi(R) is the 
volume of the window, and d2(r; R) is the volume common to two spherical windows of radius 
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R whose centers are separated by a distance r divided by t>i(i£). We will call a2{r; R) the 
scaled intersection volume. The lower bound (}3"3"j) provides another integral condition on the 
pair correlation function. 

For large R, it has been proved that c 2 (R) cannot grow more slowly than 7i? d_1 , where 
7 is a positive constant jBec87j . This implies that the Yamada lower bound in (}3"3"j) is always 
satisfied for for sufficiently large R for any d > 2. In fact, we have not been able to construct 
any examples of a pair correlation function (?2( r ) at some number density p that satisfy the 
two nonnegativity conditions (fHUJ) and (pTTj) and simultaneously violate the Yamada condition 
for any R and any d > 2. Thus, it appears that the Yamada condition is most relevant 
in one dimension, especially in those cases when cr 2 (R) is bounded. We note that point 
processes (translationally invariant or not) for which o~ 2 (R) ~ R d_1 for large R are examples 
of hyperuniform point patterns jTS03j . This classification includes all lattices as well as 
aperiodic point patterns. Hyperuniformity implies that the structure factor S'(k) has the 
following small k behavior: 

lim S(k) = 0. (34) 



The scaled intersection volume a^r; R) appearing in (J33)) will play a prominent role in 
this paper. It has support in the interval [0,2R), the range [0, 1] and the following integral 
representation: 

pcos- 1 {r/{2R)) 

a 2 (r; R) = c(d) sm d (6) d6, (35) 

Jo 

where c(d) is the <i-dimensional constant given by 

2r(l + d/2) , . 

cM = ^r ( (d + i)/2) - < 36 » 

Following the analysis given by Torquato and Stillinger |TS03] for low dimensions, we obtain 
the following new series representation of the scaled intersection volume a2{r; R) for r < 2R 
and for any d: 



«„ R) - 1 - + c M (37) 

where x = r/(2R). This is easily proved with the help of Maple. For even dimensions, relation 
(|3"Tj) is an infinite series, but for odd dimensions, the series truncates such that «2(r; R) is 
a univariate polynomial of degree d. Except for the first term of unity, all of the terms in 
(}3~Tj) involve only odd powers of x. Figure ^ shows graphs of the scaled intersection volume 
ct2(r; R) as a function of r for the first five space dimensions. For any dimension, a(r; R) is a 
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monotonically decreasing function of r. At a fixed value of r in the interval (0, 2R), 0:2(7"; R) 
is a monotonically decreasing function of the dimension d. For large d, we will subsequently 
make use of the asymptotic result 

<*2&R)~{*) [I) ^ (38) 



7T 




r/(2R) 



Figure 1: The scaled intersection volume «2( r ;-R) for spherical windows of radius R as a 
function of r for the first five space dimensions. The uppermost curve is for d = 1 and 
lowermost curve is for cf = 5. 

Before closing this section, it is useful to note that there has been some recent work that 
demonstrates the existence of point processes for a specific p and g 2 provided that p and g 2 
meet certain restrictions. For example, Ambartzumian and Sukiasian proved the existence 
of point processes that come from Gibbs measures for a special #2 for sufficiently small p 
[AS91 . Determinantal point processes have been considered by Soshnikov SosOOj and Costin 
and Lebowitz CL04 . Costin and Lebowitz have also studied certain one- dimensional renewal 
point processes |CL04j . Stillinger and Torquato |ST04j discussed the possible existence of 
a general interparticle pair potential (associated with a Gibbs measure) for a given p and 
g 2 using a cluster expansion procedure but did not address the issue of convergence of this 
expansion. Koralov |Kor05j indeed proves the existence of a pair potential on a lattice (with 
the restriction of single occupancy per lattice site) for which p is the density and g 2 is the 
pair correlation function for sufficiently small p and g 2 . There is no reason to believe that 
Koralov's proof is not directly extendable to the case of a point process corresponding to 
a sphere packing in R d , where the nonoverlap condition is the analog of single occupancy 
on the lattice. Thus, we expect that one can prove the existence of a pair potential in 
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corresponding to a sphere packing for a given p and g 2 provided that p and g 2 sufficiently 
small. 

4 Disordered Packings in High Dimensions and the Decor- 
relation Principle 

In this section, we examine the asymptotic behavior of certain disordered packings in high 
dimensions and show that unconstrained spatial correlations vanish asymptotically, yielding 
a decorrelation principle. We define a disordered packing in M. d to be one in which the pair 
correlation function #2(1") decays to its long-range value of unity faster than | r | rf e for some 
e > 0. The decorrelation principle as well as a number of other results (that will be discussed 
in Section 5) motivate us to propose a conjecture in Section 5 that describes the circumstances 
in which the two standard nonnegativity conditions (jrSOj) and (f?TTJ) are necessary and sufficient 
to ensure the existence of a disordered sphere packing. 

4.1 Example 1: Disordered Sequential Packings 

First we show that there exists a disordered sphere packing that realizes the greedy lower 
bound © (0 = l/2 d ) for all d. Then we study the asymptotic properties of the n-particle 
correlation functions in the large dimension limit. 

The disordered packing that achieves the greedy lower bound is a special case of a general- 
ization of the so-called random sequential addition (RSA) process |Tor02j . This generalization, 
which we introduce here, is a subset of the Poisson point process. Specifically, the centers of 
"test" spheres of unit diameter arrive continually throughout during time t > according 
to a translationally invariant Poisson process of density per unit time 77, i.e., rj is the number 
of points per unit volume and time. Therefore, the expected number of centers in a region 
of volume Q during time t is rjflt and the probability that this region is empty of centers is 
exp(— rjQt). However, this Poisson distribution of test spheres is not a packing because the 
spheres can overlap. To create a packing from this point process, one must remove test spheres 
such that no sphere center can lie within a spherical region of unit radius from any sphere 
center. Without loss of generality, we will set 77 = 1. 

There is a variety of ways of achieving this "thinning" process such that the subset of 
points correspond to a sphere packing. One obvious rule is to retain a test sphere at time t only 
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if it does not overlap a sphere that was successfully added to the packing at an earlier time. 
This criterion defines the well-known RSA process in M. d |Tor02| . and is clearly a statistically 



homogeneous and isotropic sphere packing in M. d with a time-dependent density <f)(t). In the 
limit t — > oo, the RSA process corresponds to a saturated packing with a maximal or saturation 
density 4> s (oc) = lim^ oo 0(t). In one dimension, the RSA process is commonly known as the 
"car parking problem", which Rehyi showed has a saturation density </> s (oo) = 0.7476... 
|Reh63j . For 2 < al < oo, an exact determination of </> s (oo) is not possible, but estimates for 
it have been obtained via computer experiments for low dimensions Tor02j. However, as we 
will discuss below, the standard RSA process at small times (or, equivalently, small densities) 
can be analyzed exactly. 

Another thinning criterion retains a test sphere centered at position r at time t if no other 
test sphere is within a unit radial distance from r for the time interval nt prior to t, where k 
is a positive constant in the interval [0, 1]. This packing is a subset of the RSA packing, and 
therefore we refer to it as the generalized RSA process. Note that when k — 0, the standard 
RSA process is recovered, and when k = 1, a relatively unknown model due to Matern [Mat86j 
is recovered. The latter is amenable to exact analysis. 

The time-dependent density <f>(t) in the case of the generalized RSA process with k = 1 
is easily obtained. (Note that for any < k < 1, the generalized RSA process is always an 
unsaturated packing.) In this packing, a test sphere at time t is accepted only if does not 
overlap an existing sphere in the packing as well as any previously rejected test sphere (which 
we will call "ghost" spheres.) An overlap cannot occur if a test sphere is outside a unit radius 
of any successfully added sphere or ghost sphere. Because of the underlying Poisson process, 
the probability that a trial sphere is retained at time t is given by exp(— vi(l)t), where t>i(l) is 
the volume of a sphere of unit radius having the same center as the retained sphere of radius 
1/2. Therefore, the expected number density p(t) and packing density 4>(t) at any time t are 
respectively given by 

p(t)= [ expi-v^l^alf = -^—il-exp(-v 1 (l)t)} (39) 
Jo v i\ L ) 

and 

4>{t) = p(tM(l/2) = 1[1 - exp(- % (l)t)]. (40) 

We see that <f>(t) is a monotonically increasing function of t. This result was first given by 
Matern using a different approach and he also gave a formal expression for the time-dependent 
radial distribution function g2{f]t) (see Section^. Here we present an explicit expression for 
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g 2 (r; t) at time t for any dimension d: 

G(r- 1) 



2^[f3 2 (r;l)-l]^(t) 



2 d <p(t) 



l-exp[-2 d /3 2 (r;l)£] 
fc{r) 1) 



Here 



Q(x) 



0, x < 0, 

1, x > 0, 



is the unit step function and 



(5 2 {r-R) = 2-a 2 (r-R) 



(41) 

(42) 
(43) 



is the union volume of two spheres of radius R (whose centers are separated by a distance r) 
divided by the volume of a sphere of radius R and a2{r; R) is the scaled intersection volume 
of two such spheres given by ([35)1 . Our approach for obtaining (|41|) is different than Matern's 
and details are given elsewhere [TS06 . 



It is useful to note that at small times or, equivalently, low densities, formula (}4T)|) yields 
the asymptotic expansion <f)(t) = v\(l)t/2 d — v f(l)t 2 /2 d+l + 0(t 3 ), which when inverted yields 
t = 2 d <p/v 1 (i) + 2 d - 1 4> 2 + 0(4> 3 ). Substitution of this last result into JUJ) gives 



g 2 (r ] <f ) ) = Q(r-l) + 0(<f ) 3 ), 
which implies that g 2 {T; 4>) tends to the unit step function 0(r — 1) as 



(44) 

for any finite d. 
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Figure 2: Radial distribution function for the first five space dimensions at the maximum 
density (p = l/2 d for the generalized RSA model with k — 1. 
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In the limit t — > oo, the maximum density is given by 

0(oo)= hm0(t) = l (45) 

and 

x , / x 20(r-l) 6(r-l) 
j2 (r; oo) S hm j2 (r; f) = - J-^ = , _ ^./^ (46) 

We see that the greedy lower-bound limit on the density is achieved in the infinite-time limit 
for this sequential but unsaturated packing. This is the first time that such an observation 
has been made. Obviously, for any < k < 1, the maximum (infinite-time) density of the 
generalized RSA packing is bounded from below by l/2 d (the maximum density for k — 1). 
Note also that because /^(r; 1) is equal to 2 for r > 2, gz(r\ oo) = 1 for r > 2, i.e., spatial 
correlations vanish identically for all pair distances except those in the small interval [0,2). 
Even the positive correlations exhibited for 1 < r < 2 are rather weak and decrease with 
increasing dimension. The function g 2 (r\ oo) achieves its largest value at r = 1 + in any 
dimension and for d — 1, g2(l + ; oo) = 4/3. The radial distribution function g2(r;oo) is 
plotted in Fig. |21 for the first five space dimensions. Using the asymptotic result (|38|) and 
relation (J46j) . it is seen that for large d, 



^ oo) : , (47) 

i_ ±\ i° 



2tt J \AJ d 1 / 2 

and thus g^ir] oo) tends to the unit step function 9(r — 1) exponentially fast as d — > oo 
because the scaled intersection volume Q2(l; 1) vanishes exponentially fast. 

The higher-order correlation functions for this model have not been given previously. In 
another work [TS06J, we use an approach different from the one used by Matern to obtain not 
only g 2 but an explicit formula for the general n-particle correlation function g n , defined by 
(|21j) . for any time t and n and for arbitrary dimension d. To our knowledge, this represents the 
first exactly solvable disordered sphere packing model for any d. These details are somewhat 
tangential to the present work and for our purposes it suffices to state the final result in the 
limit t — > oo for n > 2: 



i<j 

Pn( r l, ■ ■ ■ j r n] 1) 



<)„ ( r i . . . . , r n ; oo) = — ^ | g n -i(Qu oo) (48) 

i=i 



where the sum is over all the n distinguishable ways of choosing n—1 positions from n positions 
ri, . . . r n and the arguments of g n -\ are the associated n — 1 positions, which we denote by Qi, 
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and #i = l. Moreover, /3 n (r 1; . . . , r n ; R) is the union volume of n congruent spheres of radius 
R, whose centers are located at r 1; . . . , r n , where = \rj — r$| for all 1 < i < j < n, divided 
by the volume of a sphere of radius R. 

Lemma 4.1 In the limit d — *> oo, the n-particle correlation function g n {^i, ■ ■ ■ > r nj °°) ~ 1 
uniformly in (r 1; . . . , r n ) e M. d such that r^ > 1 for all 1 < % < j < n. If < 1 for any pair 
of points Ti and Tj, then g n ( r i, ■ ■ ■ > r n] oo) = 0. 

Proof. The second part of the Lemma is the trivial requirement for a packing. Whenever 
fij > 1 for all 1 < i < j < n, it is clear from ()48)1 that we have the following upper and lower 
bounds on the n-particle correlation function: 

£ < * < («) 

where denotes the largest possible value of g n -i- The scaled union volume j3 n of n spheres 
obeys the bounds 

n - ^2 a ^ r iv l ) ^ Pn < n, (50) 

i<j 

but since the scaled intersection volume of two spheres ct2{r; 1) attains its maximum value for 
r > 1 when r = 1, we also have 

n(n — 1) 



n — 



-a 2 (l; 1) < Pn < n. (51) 



2 

Use of this inequality and the recursive relation (}4*5j) yields the bounds 

i<g n < ? tt-^t^ ■ ( 52 ) 

Using the asymptotic result (|3~%j). we see that the upper bound tends to the lower bound for 
any given n as d — » oo, which proves the Lemma. 

In summary, the Lemma enables us to conclude that in the limit d — > oo and for = l/2 d 

n 

g n (r 12 , . . . ,r ln ; oo) ~ JJflijfoi; oo), (53) 

i<j 

where 

<7 2 (r;oo)~e(r-l). (54) 

Importantly, we see that the asymptotic behavior of g 2 in the low-density limit 
for any d [cf. f!44|) ] is the same as the high- dimensional limit d — > oo [cf. (I54J) ]. i.e., spatial 
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correlations, which exist for positive densities at fixed d, vanish for pair distances beyond 
the hard-core diameter. Note also that g n for n > 3 asymptotically factorizes into products 
involving only the pair correlation function g 2 . Is the similarity between the low-density 
and high- dimensional limits for this model of a disordered packing a general characteristic of 
disordered packings? In what follows, we discuss another disordered packing that has this 
attribute and subsequently formulate what we refer to as a "decorrelation principle." 



4.2 Example 2: The Classic Gibbsian Hard-Sphere Packing 

The statistical mechanics of the classic Gibbsian hard-sphere packing is well established (see 
[Tor02 j and references therein). The purpose of this subsection is simply to collect some results 
that motivate the decorrelation principle. Let <&n(t n ) be the iV-body interaction potential 
for a finite but large number of particles with configuration = {1*1, r 2 , . . . , r^r} in a volume 
V in at absolute temperature T. A large collection of such systems in which N, V and 
T are fixed but in which the particle configurations are otherwise free to vary is called the 
Gibbs canonical ensemble. Our interest is in the thermodynamic limit, i.e., the distinguished 
limit in which N — > 00 and V — > 00 such that the number density p = N/V exists. For a 
Gibbs canonical ensemble, when the n-particle densities p n (defined in Section EJ) exist, they 
are entirely determined by the interaction potential ^n(t n ). For a hard-sphere packing, the 
interaction potential is given by a sum of pairwise terms such that 

N 

^(^) = ^n 2 (|r,- rj |). (55) 
where Ui(r) is the pair potential defined by 

( \ \ + °° ' r<1 > fKR\ 

u 2 (r) = < (56) 
[0, r > 1, 

Thus, the particles do not interact for interparticle separation distances greater than or equal 
to unity but experience an infinite repulsive force for distances less than unity. The hard 
spheres have kinetic energy, and therefore a temperature, but the temperature enters in a 
trivial way because the configurational correlations between the spheres are independent of the 
temperature. We call this the classic equilibrium sphere packing, which is both translationally 
and rotationally invariant. 

In one dimension, the n-particle densities p n for such packings are known exactly in 
the thermodynamic limit. The density <fi lies in the interval [0, 1] but this one-dimensional 
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Figure 3: The radial distribution function for the classic three-dimensional equilibrium packing 
at <p = 0.49 as obtained from molecular dynamics computer simulations. The graph is adapted 
from Figure 3.15 of Torquato [Tor02j. 

packing is devoid of a discontinuous (first-order) transition from a disordered (liquid) phase to 
an ordered (solid) phase. Although a rigorous proof for the existence of a liquid-to-solid phase 
transition in two or three dimensions is not yet available, there is overwhelming numerical 
evidence (as obtained from computer simulations) that such a transformation takes place at 
sufficiently high densities. The maximal densities for equilibrium sphere packings in two and 
three dimensions are max = 7t/a/T2 and max = 7r/\/l8, respectively, i.e., they correspond to 
the density of the densest sphere packing in the respective dimension. 

Figure H3 shows the three-dimensional radial distribution function as obtained from com- 
puter simulations for a density = 0.49, which is near the maximum value for the stable 
disordered branch. It is seen that the packing exhibits short-range order [i.e., g2{r) has both 
positive and negative correlations for small r], but (72 (r) decays to its long-range value expo- 
nentially fast after several diameters. By contrast, in the limit d — *■ 00, it has been shown that 
the "pressure" |Rue99j of an equilibrium packing is exactly given by the first two terms of its 
asymptotic low-density expansion for some positive density interval [0, 0o] [WRF87, FP99 . 
(Roughly speaking, the pressure is the average force per unit area acting on an "imaginary 
planar wall" in the packing due collisions between the spheres and the wall.) Frisch and Percus 
[FP99J have established, albeit not rigorously, that O — l/2 d . This result for the pressure 
implies that the leading-order term of the low-density expansion of the radial distribution 
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function in arbitrary dimension |Tor02j 

g 2 ( r ) = e(r - 1) [l + 2%(r; 1)0 + 0(0 2 )] (57) 

becomes asymptotically exact in the limit d —>■ oo in the same density interval. The presence 
of the unit step function 0(r — 1) in relation (|57j) means that the scaled intersection volume 
a 2 (r; 1) need only be considered for values of r in the interval [1,2]. Since a 2 (r; 1) is largest 
when r = 1 for 1 < r < 2 and a 2 (l; 1) has the asymptotic behavior (|38jh the product 
2 d a 2 (l; 1)0 vanishes no more slowly than (6/7r) 1 / 2 /[(4/3) d / 2 rf 1 / 2 ] in the limit d -> oo for < 
< l/2 d , and therefore g 2 (r) tends to @(r — 1) exponentially fast. In summary, we see again 
that spatial correlations that exist in low dimensions for r > 1 completely vanish in the limit 
d — > oo. Moreover, this is yet another disordered packing model in which the high-dimensional 
asymptotic behavior corresponds to the low-density asymptotic behavior. 

The corresponding n-particle correlation function g n , defined by (|2*T|) . in the low-density 
limit Sal58 is given by 



g n {r 12 ,...,r ln ) = JJ g 2 (r ij ) [l + 2 d a n {r 12 , . . . , r ln ; 1)0 + 0(0 2 )] (58) 

i<j 

where a n (ri 2 , . . . , i"i„; R) is the intersection volume of n congruent spheres of radius R (whose 
centers are located at ri, . . . , r n , where r»j = rj — r.j for all 1 < % < j < n) divided by the 
volume of a sphere of radius R. The scaled intersection volume a n (r 12 , . . . , r ln ; R)/n has the 
range [0,1]. Now since a^r^-, 1) > a„(r 12 , . . . , r ln ; 1) for any pair distance = [ r^- 1 such 
that 1 < i < j < n, then it follows from the analysis above that in the limit d — > oo for 
< < l/2 d 

g 2 (r) ~9(r-l) (59) 

and 

n 

g n (r 12 , . . . ,n n ) ~ JJ^^ij)- (60) 

Again, as in the generalized RSA example with k — 1, g n factorizes into products involving 
only g 2 s in the limit d — *■ oo. Moreover, we should also note that the standard RSA process 
(generalized RSA process with k — 0) has precisely the same asymptotic low-density behavior 
as the standard Gibbs hard-sphere model [Tor02 . More precisely, these two models share 
the same low-density expansions of the g n through terms of order and therefore the same 
asymptotic expressions (jo^j) and (|5Uj) . 
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4.3 Decorrelation Principle 



The previous two examples illustrate two important and related asymptotic properties that 
are expected to apply to all disordered packings: 

1. the high-dimensional asymptotic behavior of g 2 is the same as the asymptotic behavior 
in the low-density limit for any finite d, i.e., unconstrained spatial correlations, which 
exist for positive densities at fixed d, vanish asymptotically for pair distances beyond 
the hard-core diameter in the high- dimensional limit; 

2. and g n for n > 3 asymptotically can be inferred from a knowledge of only the pair 
correlation function g 2 and number density p. 

What is the explanation for these two related asymptotic properties? Because we know from 
the Kabatiansky and Levenshtein (1979) asymptotic upper bound on the maximal density that 
4> must go to zero at least as fast as 2~ ' 599 for large d, unconstrained spatial correlations 
between spheres must vanish, i.e., statistical independence is established. (An example of 
constrained spatial correlations is described below.) Such a decorrelation means that the g n 
for n > 3 are determined entirely from a knowledge of the decorrelated pair correlation function 
#2- In the specific examples that we considered, the g n factorize into products involving only 
<72 ! s, but there may be other decompositions. For example, the g n for n > 3 can be functionals 
that only involve p and g 2 . We will call the two asymptotic properties the decorrelation 
principle for disordered packings. This principle as well as other results described in Section 
5 leads us to a conjecture concerning the existence of disordered sphere packings in high 
dimensions, which we state in Section 5.1. 

An example of constrained spatial correlations that would not vanish asymptotically is 
illustrated in Fig. 01 where we show the pair correlation function g 2 {r) for a three-dimensional 
sphere packing near the so-called maximally random jammed state |TTD00[ ITor02j . A special 
feature of <72( r ) for a maximally random jammed packing is a delta-function contribution 
at r = 1, which reflects the fact that the average kissing number (i.e., average number of 
contacting particles per particle) is effectively six for this collectively jammed packing, meaning 
that the packing is isostatic |DST05j . A positive average kissing number is required if the 
packing is constrained to be jammed and in M. d this means that the average kissing number 
is 2d for either collective or strict jamming [DST05 . Isostatic packings are jammed packings 
with the minimum number of contacts for a particular jamming category. According to the 
Decorrelation Principle, as d tends to infinity, g 2 for a maximally random jammed packing 
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Figure 4: The radial distribution function of a three-dimensional packing of spheres near 
the maximally random jammed state |TTD00t ITor02j at a density <fi = 0.64 as obtained 



from computer simulations |TS02j . The delta function contribution at r = 1 (of course, not 
explicitly shown) corresponds to an average kissing number of about six. 

would retain this delta-function contribution but the unconstrained spatial correlations beyond 
r = 1 would vanish. Of course, the manner in which the g% shown in Fig. ^approaches the 
asymptotic limit of a step function 0(r — 1) plus a delta-function contribution at r = 1 is 
crucial. We note that maximally random jammed packings contain about 2-3% rattler spheres, 
i.e., spheres trapped in a cage of jammed neighbors but free to move within the cage. 



5 New Approach to Lower Bounds 

The salient ideas behind our new approach to the derivation of lower bounds on max were 
actually laid out in our earlier work |TS02j . The main objective of that work was to study 
sphere packings in three dimensions in which long-range order was suppressed and short- 
range order was controlled (i.e., disordered sphere packings in 9ft 3 ) using so-called (72-invariant 
processes. A g2-invariant process is one in which a given nonnegative pair correlation ^(r) 
function remains invariant for all r over the range of densities 

O<0<0*. (61) 

The terminal density is the maximum achievable density for the gi~ invariant process subject 
to satisfaction of the structure factor S*(k) inequality (J3TJ) . A five-parameter test family of 
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g 2 s had been considered, which incorporated the known features of core exclusion, contact 
pairs, and damped oscillatory short-range order beyond contact. The problem of finding the 
maximal packing fraction 0* was posed as an optimization problem: maximize over the set 
of parameters subject to the constraints (fHUJ) and (pTTj) . We noted in passing that when the 
damped-oscillatory contribution to g 2 was set equal to zero, the optimization problem could 
be solved analytically for all space dimensions, leading to a terminal density 0* = (d + 2)/2 d+1 . 
Under the assumption that such a g 2 was a realizable packing, we also observed that this 0* 
was a lower bound on the maximal density for any sphere packing [i.e., max > (d + 2)/2 d+1 ] 
because the terminal density would have been higher by including the damped-oscillatory 
contribution to g 2 . This conjectural lower bound was noted to provide linear improvement 
over Minkowski's lower bound, but we were not aware of Ball's similar lower bound [Bal92j 
at the time. Since our original 2002 paper, we also learned about other necessary conditions 
for the realizability of a point process for a given number density p and g 2 , such as Yamada's 
condition (|H2*|) . In any event, our brief remarks about lower bounds on sphere packings were 
not intended to be mathematically rigorous. 

It is our intent here to make our optimization methodology to obtain lower bounds on max 
more mathematically precise, especially in light of recent developments and the considerations 
of the previous two sections. We then apply the optimization procedure to provide alternative 
derivations of previous lower bounds as well as a new bound. 

We will consider those "test" g^r)^ that are distributions on ~R d depending only on the 
radial distance r such that h(r) = g 2 (r) — 1. For any test g 2 (r), we want to maximize the 
corresponding density satisfying the following three conditions: 

(i) S^O") > for all r, 

(ii) supp(# 2 ) C {r : r > 1}, 
(ih) 

S(k) = l+p (2ix)i f°° r^lg^r) - 1] J(d/2) ,~ff dr > for all k. 
Jo (kry 

We will call the maximum density the terminal density and denote it by 0*. 
Remarks: 

1. The conditions (i)-(iii) are just recapitulations of (|2*5Jl . (fHUj) . and (pTTj) for this class of test 
functions. We will call condition (ii) the hard-core constraint. 

2. When there exist sphere packings with g 2 satisfying conditions (i)-(iii) for in the interval 
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[0, </>*] , then we have the lower bound on the maximal density given by 

0max > </>*• (62) 

The best lower bound would be obtained if one could probe the entire class of test functions. 
In practice, we will consider here only a small subset of test functions and in particular those 
that are amenable to exact asymptotic analysis. In some instances, we will associate with the 
terminal density 0* an optimized average kissing number Z*. Thus, whenever inequality (f6*2*j) 
applies, the maximal kissing number Z max is bounded from below by Z*, i.e., 

^max > Z*- (63) 

In the next subsection, we put forth a conjecture that states when the conditions (i)-(iii) are 
necessary and sufficient for the existence of disordered sphere packings. 

3. Remarkably, the optimization problem defined above is identical to one formulated by Cohn 
|Coh02j . In particular, it is the dual of the primal infinite-dimensional linear program that 
Cohn employed with Elkies |CE03j to obtain upper bounds on the maximal packing density. 
One need only replace S(k) with g — c5(k), where c plays the role of number density, g is a 
tempered distribution and g is its Fourier transform in Cohn's notation. Thus, even if there 
does not exist a sphere packing with g 2 satisfying conditions (i)-(iii), our formulation has 
implications for upper bounds on max , which we discuss in Section |B] For finite-dimensional 
linear programs (and many infinite-dimensional ones) there is no "duality gap," i.e., the optima 
of the primal and dual programs are equal. However, in this infinite-dimensional setting, it 
is not clear how to prove that there is no duality gap |Coh02j . Therefore, it is rigorously 
true that the terminal density 0* can never exceed the Cohn-Elkies upper bound, which is a 
desirable feature of our formulation, otherwise the terminal density could never correspond to 
a rigorous lower bound. 

We will show that for the test radial distribution functions considered in this paper, the 
Yamada condition (}3*2*|) is only relevant in one dimension, and even then in just some cases. 
A remark about the Yamada condition for sphere packings is in order here. In earlier work 
[TS03J, we observed that for any sphere packing of congruent spheres, the number variance 
for a spherical window of radius R defined by (fBlIj) obeys the lower bound 

(64) 

for any R. This is a tight bound for sufficiently small R and is exact for R < 1/2. However, 
we note here that provided that 2 d <ftR d < 1, the Yamada lower bound (}3*2*j) and lower bound 



a 2 (R) > 2 d (t)R d 



1 - 2 d cj)R d 
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(jSljl are identical. Thus, the Yamada lower bound for any sphere packing only needs to be 
checked for R > Rq, where 

fl » = ^Af < 65 > 

5.1 Existence of Disordered Packings in High Dimensions 

We have seen that a necessary condition for the existence of a translationally invariant point 
process with a specified positive p and nonnegative g 2 is that S(k) is nonnegative [cf. (|3ip] . 
In other words, given p and g 2 , it does not mean that there are some higher-order functions 
<73,<?4, ... for which these one- and two-point correlation functions hold. The function g 2 
specifies how frequently pair distances of a given length occur statistically in M d . The third- 
order function g 3 reveals how these pair separations are linked into triangles. This additional 
information generally cannot be inferred from the knowledge of p and g 2 alone, however. 
The fourth-order function g± controls the assembly of triangles into tetrahedra (and is the 
lowest-order correlation function that is sensitive to chirality) but g^ cannot be determined 
by only knowing p, g 2 and g%. In general, g n for any n > 3 is not completely determined from 
a knowledge of the lower-order correlations functions alone. This is to be contrasted with 
general stochastic processes in which nonnegativity of first- and second-order statistics (mean 
and autocovariance) are necessary and sufficient to establish existence because one can always 
find a Gaussian process with such given first- and second-order statistics. For a Gaussian 
process, first- and second-order statistics determine all of the high-order statistics. 

There are a number of results that suggest it is reasonable to conclude that the gener- 
ally necessary nonnegativity conditions for the existence of a disordered sphere packing be- 
come necessary and sufficient for sufficiently large d. First, the decorrelation principle of the 
previous section states that unconstrained correlations in disordered sphere packings vanish 
asymptotically in high dimensions and that the g n for any n > 3 can be inferred entirely from 
a knowledge of p and g 2 . Second, as we noted in Section^ the necessary Yamada condition 
appears to only have relevance in very low dimensions. Third, we will demonstrate below that 
other new necessary conditions also seem to be germane only in very low dimensions. Fourth, 
we will describe numerical constructions of configurations of disordered sphere packings on the 
torus corresponding to certain test radial distributions functions in low dimensions for den- 
sities up to the terminal density. Finally, we will show that certain test radial distributions 
functions recover the asymptotic forms of known rigorous bounds. In light of these results, 
we propose the following conjecture: 
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Conjecture 5.1 For sufficiently large d, a hard-core nonnegative tempered distribution <?2( r ) 
that satisfies <?2( r ) = 1 + 0(|r| -d-e ) for some e > is a pair correlation function of a 
translationally invariant disordered sphere packing in M d at number density p if and only 
if S(\t) = 1 + ph(k) > 0. The maximum achievable density is the terminal density 0*. 

Remarks: 

1. A weaker form of this conjecture would replace the phrase "for sufficiently large d" with 
"in the limit d — > oo." 

2. Employing the aforementioned optimization procedure with a certain test function g 2 and 
this conjecture, we obtain in what follows conjectural lower bounds that yield the long-sought 
asymptotic exponential improvement on Minkowski's bound. Before obtaining this result, we 
first apply the procedure to two simpler test functions that we examined in the past. 



5.2 Step Function 

The simplest possible choice for a radial distribution function corresponding to a disordered 
packing is the following unit step function: 

02 (r) = 0(r - 1). (66) 

This states that all pair distances beyond the hard-core diameter are equally probable, i.e., 
spatial correlations vanish identically. The corresponding structure factor [cf. condition (iii)] 
for this test function in any dimension d is given by [TS02 

m _ 1 _ ^±A Mk) , m 

where v = d/2. Since there are no parameters to be optimized here, the terminal density <fi* is 
readily obtained by determining the highest density for which the condition (|5Tj) is satisfied, 
yielding 

^ = ( 68 ) 



Now we show that the Yamada condition (JHzj) is satisfied in any dimension for < (f> < 
2~ d . Consider the more general class of radial distribution functions: 

< g 2 {r) < 1 for r > 1. (69) 
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The test function (jbbj) belongs to this class. Note that for any dimension, the scaled intersec- 
tion volume given by (|37jl obeys the inequality 

a 2 (r; R)<l-S- for < r < 2R, (70) 
2R 

where the equality applies when d — 1. For (?2( r ) satisfying (joTJj) . relation (jH^j) and inequality 
(f7U|) yield the following lower bound for any <i: 



a 2 (R) > 2 l 



l + d2 d (j) / r d ~ l h{r) ' 



o 



1 - 

2R 



dr 



(71) 



At = l/2 d , the lower bound (f7Tj) for the test function (JoT)J) is given by 

ff2(fl) £ 2(5TT) fld " (72) 

and because i? — 1 [cf- we only need to consider i? > 1. In particular, the right side of 

this inequality is smallest at R — 1 so that 

^ £ 2I5TT)' < 73 ' 

Since u 2 {R) > 1/4 for d > 1, Yamada's condition is satisfied for all R for the step function 
fU) at (j) = 1 /2 d as well as all < l/2 d . 

We already established in Section 0] that there exist sphere packings that asymptotically 
have radial distribution functions given by the simple unit step function (}66|) for < 2~ d . 
Nonetheless, invoking Coni ecture 15 . 1 1 and terminal density specified by (J68|) implies the asymp- 
totic lower bound on the maximal density is given by 

0max > ^, (74) 

which provides an alternate derivation of the elementary bound (J7|). 

Using numerical simulations with a finite but large number of spheres on the torus, we 
have been able to construct particle configurations in which the radial distribution function 
(sampled at discretized pair distances) is given by the test function in one, two and 
three dimensions for densities up to the terminal density |(JTS031 IUST06] . The existence 
of such a discrete approximation to ()66|) of course is not conclusive proof of the existence 
of such packings in low dimensions, but they are suggestive that the standard nonnegativity 
conditions may be sufficient to establish existence in this case for densities up to 0*. 
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5.3 Step Plus Delta Function 



An important feature of any dense packing is that the particles form contacts with one another. 
Ideally, one would like to enforce strict jamming (see Introduction). The probability that a 
pair of particles form such contacts at the pair distance r = 1 for the test function (jfifij) is 
strictly zero. Accordingly, let us now consider the test radial distribution function given by 
the previous test function plus a delta-function contribution as follows: 

^( r ) = e(r-l) + -^-5(r-l). (75) 

Here si(r) is the surface area of a ci-dimensional sphere of radius r given by (ffifj) and Z is a 
parameter, which is the average kissing number. Because we allow for interparticle contacts 
via the second term in l)75j) . the terminal density is expected to be greater than 2~ d , which 
will indeed be the case. The corresponding structure factor [cf. (hi)] for this test function in 
any dimension d is given by [TS02 



ff( *) = i -^±A m + ^Aj^ {k) , (76) 



where v = d/2. The structure factor for small k can be expanded in a MacLaurin series as 
follows: 



S(k) = 1 + (Z- 2 d <t>) + 



2^ Z 



k 2 + 0(k 4 ). (77) 



l + d/2 2d 

The last term changes sign if Z increases past 2 d (pd/(d+ 2). At this crossover point, 

2<m 

S(k) = l- JT ^^ + 0(k 4 ) (78) 

Under the constraint that the minimum of S(k) occurs at k — 0, then we have the exact 
results 

d + 2 d 

We see that at the terminal density, the average kissing number Z* = d/2, which does not 
even meet the local jamming criterion described in the Introduction. 

The Yamada condition (J32)) is violated only for d — 1 for the test function (J75|) at the 
terminal density specified by (fT^j). It is easy to verify directly that the Yamada condition 
becomes less restrictive as the dimension increases from d = 2. Interestingly, we have also 
shown via numerical simulations that there exist sphere packings possessing radial distribution 
functions given by the test function (fTH|) (in the discrete approximation) in two and three 
dimensions for densities up to the terminal density [OST06J. This is suggestive that Conjecture 
15.11 for this test function may in fact be stronger than is required. 
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In the high-dimensional limit, we invoke Conjecture 15.11 and the terminal density given 
by (|79|h yielding the conjectural lower bound 

d + 2 , . 

0max > -^p (80) 

This lower bound provides the same type of linear improvement over Minkowski's lower bound 
as does Ball's lower bound jBal92j. 



5.4 Step Plus Delta Function with a Gap 

The previous test function (|75j) provided an optimal average kissing number Z* = d/2 that did 
not even meet the local jamming criterion. Experience with disordered jammed packings in low 
dimensions reveals that the kissing number as well as the density can be substantially increased 
if there is there is a low probability of finding noncontacting particles from a typical particle at 
radial distances just larger than the nearest-neighbor distance. This small-distance negative 
correlation is clearly manifested in the graph of (r) for the three-dimensional maximally 
random jammed packing (Figure HJ for values of r approximately between 1.1 and 1.5. We 
would like to idealize this small- distance negative correlation in such a way that it is amenable 
to exact asymptotic analysis. Accordingly, we consider a test radial distribution function that 
is similar to the previous one [cf. (J75|) ] but one in which there is a gap between the location 
of the unit step function and the delta function at finite d, i.e., 

g 2 (r) = 0(r - a) + -^r-S(r - 1). (81) 

The expression contains two adjustable parameters, a > 1 and Z, which must obviously be 
constrained to be nonnegative. According to the decorrelation principle of Section 01 the 
location of the step function r = a must approach unity asymptotically, i.e., it must approach 
the previous test function (|75|). However, as we have emphasized, the manner in which the test 
function (JHTj) approaches (|73|l is crucial. Indeed, we will see that the presence of a gap between 
the unit step function and delta function will indeed lead asymptotically to substantially higher 
terminal densities. 

The structure factor is given by 

The goal now is to find the optimal values of the the adjustable nonnegative parameters Z and 
a that maximize the density <fi subject to the constraint (iii). This search in two-dimensional 
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parameter space can be reduced by imposing the further condition that a minimum of the 
structure factor occurs at k = 0. The MacLaurin expansion of expression (J%2*j) gives 



S = 1 + [Z- (2a) d <j)] + 



2<i-l„d+2, 



a 



d + 2 



Z 
2d 



k 2 + 0(k 4 



(83) 



Requiring that a zero of S(k) occurs at the origin (hyperuniformity) such that the quadratic 
coefficient is nonnegative implies the restrictions 



Z = (2a) d (f>- 1 



and 



(2a) d (f)[da 2 - (d + 2)] + d + 2 > 0. 
Combination of (J82j) and (|84j) yields the structure factor as 

Q(l\ 1 (j\ J v( k(J ) . fj\ J v-i( k ) 

s{k) = 1 - Cl{d) jkW + C2i) V^ 

where the <i-dependent coefficients ci(d) and 02(d) are given by 

Cl (d) = <Pa d 2 Zv T(l + v) 

X(l + u) 



C2(d) = [(2a) d <P-l]2»- 



d 



(84) 
(85) 

(86) 

(87) 
(88) 



Now the problem reduces to finding the optimal value of the parameter a(d) as a function of 
the space dimension d that maximizes the density <fr subject to the constraint (iii). It will be 
shown below that the optimal a is of order unity and approaches unity in the limit d — *■ 00. It 
immediately follows from (J86|) and the asymptotic properties of the Bessel functions of fixed 
order that S(k) — > 1 for k — > 00. 

In general, S(k) will possess multiple minima and thus we want to ensure that the values 
of S(k) at each of these minima are all nonnegative. To find the minima of S(k), we set its 
first derivative to zero, yielding the relation 



4 ij-J v+ i(ha) = c 2 (d)kJ u (k), 



a 



where we have used the identity 



d 

dx 



JJx) 



X" 



J u +l(x) 



X" 



(89) 



(90) 



For sufficiently small d (d < 200), the search procedure is carried out numerically using 
Maple, and is made more efficient by exploiting the fact that the minima of S(k) occur at 
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Figure 5: The optimized structure factor for d — 12 and d = 24. 

the real solutions of (jEHJ)- Figure shows the optimized structure factor for d = 12 and 
d = 24. Our numerical examination of S(k) for a wide range of d values has consistently 
shown that the first minimum for positive k is the deepest one. Although we have not proven 
this rigorously, we assume that this is a general result. 

We should note that the Yamada condition (J3"2~j) is violated only for d — 1 for the test 
function (fHTj) for the terminal density 0* and associated optimized parameters <r* and Z* 
[calculated via (|84|)]. One can again verify directly that the Yamada condition becomes less 
restrictive as the dimension increases from d = 2. However, although the test function (JHTj) 
for d = 2 with optimized parameters = 0.74803, o~* = 1.2946 and = 4.0148 satisfies the 
Yamada condition, it cannot correspond to a sphere packing because it violates local geometric 
constraints specified by o~* and Z*. Specifically, for an average kissing number of 4.0148, there 
must be particles that are in contact with at least five others. But no arrangement of the five 
exists which is consistent with the assumed pair correlation function (step plus delta function 
with a gap from 1 to 1.2946). Simple geometric considerations show that either some pairs 
of the five would be forced into the gap, or they would be restricted to fixed separations that 
would correspond to undesired delta functions beyond the gap. To our knowledge, this is 
the first example of a test radial distribution function that satisfies the two standard non- 
negativity conditions (j3Tf|) and (J3T|) and the Yamada condition (J32J) . but cannot correspond to 
a point process. Thus, there is at least one previously unarticulated necessary condition that 
has been violated in the low dimension d = 2. 

In three dimensions one obtains 0* = 0.5758254, a* = 1.246997, and Z* = 7.932582. 
The last of these requires that some nonzero fraction of the spheres must have at least 8 
contacting neighbors. We have verified that valid arrangements of both 8 and 9 contacts are 
possible, thereby avoiding the analog of the violation encountered in d — 2. As is the case 
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Table 3: Optimized parameters a*, Z* and 0*, and the ratio 



2 d+i 



d+2 



which is the relative 



improvement of the terminal density over the gapless-test-function terminal density [cf. (|79j)]. 



d 




z* 


(f>* 


2 d +V. 
d-\-2 


3 


i 




7.932582 


0.5758254 


1.842641 


4 


1 


21 2589 


13.71016 


0.4252472 


2.267985 


5 


1 


1 8RQ9Q 


21.97918 


0.3048322 


2.787037 


6 






33.53884 


0.2136444 


3.418310 


7 






49.42513 


0.1471058 


4.184343 


8 


1 

J. 




70.88348 


0.09985085 


5.112364 


24 


1 


058992 


5473.546 


8.245251e-05 


106.4095 


36 


1 


041611 


76521.15 


2.566299e-07 


928.1828 


56 


1 


028036 


4.248007e06 


1.253255e-ll 


31140.19 


60 


1 


026330 


9.179315e06 


1.674130e-12 


62262.60 


64 


1 


024823 


1.968233e07 


2.221414e-13 


124175.32 


80 


1 


020211 


3.908042e08 


6.521679e-17 


1.922982e06 


100 


1 


016421 


1.478804el0 


2.288485e-21 


5.688234e08 


125 


1 


013311 


1.246172el2 


5.610270e-27 


3.758024e09 


150 


1 


011214 


9.698081el3 


1.275632e-32 


2.319290ell 


175 


1 


009671 


7.086019el5 


2.745830e-38 


1.485866el3 


200 


1 


008510 


4.959086el7 


5.667098e-44 


9.016510el4 



with the Yamada condition ()32)1 . this additional necessary condition appears to lose relevance 
as d increases. 

The terminal density 0* and the associated optimized parameters a* and are listed 
in Table El for selected values of the space dimension between d = 3 and d = 200. Note that 
for d < 56, the terminal density lies below the density of the densest known packing. For 
d = 56, the densest arrangement is a lattice (designated by L 56 ^(M) |Neb98j ) with density 
cf) = 2.327670 x 10 -11 , which is about twice as large as 0*, as shown in the table. However, 
for d > 56, 0* can be larger than the density of the densest known arrangement. For d = 60, 
the densest known packing is again a lattice (designated by L 56j2 (M) |CS98j ) with density 
<p = 2.966747 x 10~ 13 , which is about five times smaller than 0*, as shown in the table. The 
next dimension for which data is available is d — 64, where the densest known packing is the 
iVe64 lattice |Neb98j with density = 1.326615 x 10~ 12 , which is about six times larger than 
0*. The table also reveals exponential improvement of the terminal density 0* over the one 
for the gapless case, i.e., 0* = (d + 2)/2 d+1 . The crucial question is whether such exponential 
improvement persists in the high- dimensional limit. 
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To obtain an asymptotic expression for 0* for large d, we use the fact that (2a) d <j) ^> 1, 
implying that Ci(d)/ 02(d) — > d [cf. (j57j) and (jHHj) ]. Therefore, the minima of S(k) for large 0? 
are the solutions of 

^±1^ = k -J v (k). (91) 

We see that the locations of the minima depend only on o (not on <fi). The deepest minimum 
of S(k), after the one at k = 0, is a zero and occurs at the wavenumber k = k m i n . [This 
characteristic is true in any dimension (see Figure EJ).] Therefore, S(k = k min ) = 0, 02(d) = 
Ci(d)/d, and relation (JH2J) gives the condition 

-— /\ u (k m i n ) = 1, (92) 

mm 

where 

A/7 \ Ju(k m i n <7) Jv-l(kmin) (QQ\ 

a v a 

The solution of equation (|9*2)l produces the desired optimal values of a* and 0*, where 



nun 



2^Y(l+u)afK(K 



(94) 



We find the solutions of (j91)l by linearizing each Bessel function in (j91j) around its respec- 
tive first positive zero, i.e., 

J v (x) = p l (x )(x-x )+O((x-x ) 2 ), (95) 

Ju+i(x) = f3 2 {yo)(x-y ) + 0({x-y ) 2 ), (96) 

where 

Pi(x ) = ]-[J v ^(xq) - J„ + i(x )], (97) 

foiyo) = liUvo) -JMvo)}, (98) 

and xq and yo denote the locations of the first positive zeros of J v (z) and J v+ x(z), respectively. 
Similarly, we employ the linearized form 

xJ v (x) = x Pi(x )(x — x ) + O ((x — Xq) 2 ) . (99) 

Use of these relations in (}9Tj) yields the following equation for k min : 

k mm - * - . (100) 

^a v - x Xr) - da 
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This formula provides an excellent approximation for k min . For example, for d = 200 (or 
v = 100), substitution of the exact values x = 108.8361659, y = 109.8640469 and f3 1 / f3 2 = 
1.003189733 as well as the numerical search solution a* = 1.008510 into this formula predicts 
kmin = 108.4368917. This value is to be compared to the numerical search solution of k m i n = 
108.4395. This supports the fact that the higher-order terms in the aforementioned linearized 
forms of the Bessel functions are negligibly small. Indeed, we expect that this can be rigorously 
proved, but we shall not do so here. We will assume the validity of the linearized forms in the 
asymptotics displayed below. 

For large d = 2v, we make use of the asymptotic formulas 

Xo = I , + aiZ //3 + ^ + ^ + o(_L) ) (ioi) 
yo = ^ + ^ 1/3 + 1 + ^ + ^ + v + 0( i ) ' (102) 

where the constants ax, a 2 and are explicitly given in the Appendix. For d = 200, these 
formulas predict xq = 108.8362067 and yo = 109.8640871, which are in excellent agreement 
with the exact values reported in the preceding paragraph. Using the asymptotic results given 
in the Appendix, we obtain that 

A 2 2C 2 n ( 1 \ . , 

where the constants C\ and C 2 are given explicitly in terms of the constants a x and a 2 in the 
Appendix. For d = 200, for example, this formula together with (jA-19|) provide the estimate 
(5x1(52 = 1.007122331, which is to be compared to the exact result (5 X / (5 2 = 1.006215695. 

The optimized asymptotic form for cr* is obtained by taking the derivative of both sides 
of the zero-condition (|9*2^) with respect to a and solving for a using relation ()100|) for k m i n . 
We obtain that 

v z/ 5 / 3 (is 2 
where 

qx = 0.90763589355... (105) 
is the unique positive root of xe x + e 2x — 5e x + 4 = and 

* = °.(^-2^.-10e^ + 4 + ^ + 4^.) = 

H 3e«i (2gie» - 2qx + 12 + 3e 2 ^ - 13e^) v ; 

Therefore, expression (|100j) for k min has the asymptotic form 

kmin = » + axv 1/3 + Qx + -^ + G (J-^J , (107) 
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^ = 1 + - + 7^ + ^ 7^)^ (104) 



where 



Qi 



2(?i - 1) 



-0.3860921576. 



(108) 



e« - 2 

These formulas predict a* = 1.008482538 and k min = 108.4501542, which again are in excellent 
agreement with values reported above. 

Linearizing each Bessel function appearing in (J93|) about its first positive zero and using 
the results of the Appendix yields 

f3i(XQ) (kminCT* — Xq) 



/S. u (k, n 



(109) 



< 2v 
where Pi(xq) is given by (|57|). /^(^o) = [^-2(^0) — Jv{zq)]/2, and zo is the first positive zero of 
Jnu-x- Using relations ()104|) and (|107jh and the results of the Appendix yields the asymptotic 
expansion of A„(/c m j n ): 



Di D 2 ^ ( 1 
+ — ^ + 



,2/3 ,,4/3 



5/3 



where 



D 1 



Ci(2-e 



91' 



2e"i 



£9 



C 1 [a 1 (2e^ + 6 gi e-«-7) + 3g 2 (g 1 -l)] , C 2 A 



Note also that 



3 a i i/ 1 /3 + Q 1 



1 + 



3(2 -e«i) 
£1 E 2 



Ci 



(110) 



111) 



,,!./:! f ^2/3 +C> L 



and 



where 



a 



2v 



1 + 



2q 2 q\ 2q. 



v 



2/3 



+ 



+ 



v z/ 4 / 3 V 2/5/3 

,4 A „2„ 1 /|„2 



Ei = a 2 — -E2 = — Qi^i 



— 4aj a 2 + 4a 2 



(112) 



(113) 



(114) 



and Qi is given by (JTUSJ). For d = 200, these formulas [together with the constants specified by 
(fTO. (fTUoT) . (fTTPLfrm (DT5l) and (jAJ ] predict A u (k min ) = 0.00567441932, {k min /u) u = 
3353.018128 and a1 v = 5.405924156. These values should be compared to the exact value of 
A„(k min ) = 0.00559813885, [k min jv) v = 3301.799093 and = 5.445550297. 

Thus, substituting the asymptotic relations above into the optimal expression ()94p for 
the density and invoking Conjecture 15.11 yields the conjectural lower bound 



> 



2 [ 3 - lo S2 (e)] ^-log 2 (e)ai i/V3 + ( 2qi -Q 1 ) l og2 (e) 
E 2 - 2<7 2 - D 2 /D 1 ( 1 



2Di V 7T 



u l "+ El 



V 



1/6 



,1/2 



5/6 



(115) 
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where we have used the asymptotic relation T(l + v) ~ v h '\/2itve~ v '. For d = 200, this 
asymptotic formula [together with the constants specified by (jl05|) . (jl06|) . (flllj) . HHj) , ()A-5|) 
and ()A-19|) ] predicts 0* = 5.626727001 x 10~ 44 , which is in good agreement with the numerical 
search solution of 0* = 5.667098 x 10~ 44 . Note also that the formula with k min estimated 
from (jlOOJ) yields 0* = 5.666392126 x 10 -44 , which is remarkably close to the numerical 
solution. For large d, result (|115|) yields the following dominant asymptotic formula for the 
conjectural lower bound on max : 

rfV6 _ 3.276100896 d 1 / 6 

0max _ <P* ~ 2 2 /3 £) 2$- lo ^( e )] d / 2 ~ 2 ' 7786524795 - d ' 

This putatively provides the long-sought exponential improvement on Minkowski's lower 
bound. Note that the constant Di = 0.1084878572 appearing in f)116j) is determined from 
the appropriate relation in (jlll|) using the value for qi given by ()105|) and the more refined 
estimate of C\ given by (jA-21|) . 

Substitution of the asymptotic expression ()116j) into (J%3|) and use of (j6~3"|) yields a conjec- 
tural lower bound on the maximal kissing number 

Zmax > Z* ~ 2V3e2gl ^ 1/6 2 [iog 2 (e)-i]^/2 = 40.24850787 d 1 /^ - 2213475205 -", (117) 

which applies for large d. This result is superior to the best known asymptotic lower bound on 
the maximal kissing number of 2 0,2075 |Wyn65| . Note that such a disordered packing would 
be substantially hyper static (the average kissing number is greater than Id [DST05J) and 
therefore would be appreciably different from a maximally random jammed packing TTD00, 
ITor02j . which is isostatic (see Section 4.3) and hence signficantly smaller in density. 



6 Discussion 



Our results have immediate implications for the linear programming bounds of Cohn and 
Elkies [CE03J, regardless of the validity of Conjecture 15.11 As we noted earlier, our optimiza- 
tion procedure is precisely the dual of the their primal linear programming upper bound. The 
existence of our test functions (fTS]) and (JHTj) that satisfy the conditions (i), (ii) and (iii) of 
Section 5 for densities up to the terminal density 0* narrows the duality gap [cf. ()79|) and 
([115)1 ] . In particular, inequality ()115|) provides the greatest lower bound known for the dual 
linear program. Moreover, the existence of the inequalities (180)1 and ()115j) proves that linear 
programming bounds cannot possibly match the Minkowski lower bound for any dimension d. 
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Finally, this link to the Cohn-Elkies formulation proves that the terminal density 0* can never 
exceed the Cohn-Elkies upper bound, which obviously must be true if the terminal density 
corresponds to a rigorous lower bound. 

Conjecture 15.11 concerning the existence of disordered sphere packings is plausible for a 
number of reasons: (i) the decorrelation principle of Section 4.3; (ii) the necessary Yamada 
condition appears to only have relevance in very low dimensions; (iii) other new necessary 
conditions described in Section 5.4 also seem to be germane only in very low dimensions; (iv) 
there are numerical constructions of configurations of disordered sphere packings on the torus 
corresponding to these test radial distributions functions in low dimensions for densities up to 
the terminal density [CTS03, OST06J; and (v) the test radial distributions functions and 
(|75jl recover the asymptotic forms of known rigorous bounds. Concerning the latter point, if 
Conjecture 15. II is false, it is certainly not revealed by the results produced by the test functions 
(|ni3) and ([75]) because the forms of known rigorous results, obtained using completely different 
techniques, are recovered. If Conjecture 15. H is false, why would it suddenly be revealed by the 
introduction of a gap in the test radial distribution function [cf. (|81j)] relative to (|75|1? This 
would seem to be unlikely and lends credibility to the conjecture in our view. 

Conjecture 15. l\ the particular choice (JHTj) and our optimization procedure leads to a lower 
bound on the maximal density that improves on the Minkowski bound by an exponential 
factor. Our results suggest that the densest packings in sufficiently high dimensions may be 
disordered rather than periodic, implying the existence of disordered classical ground states 
for some continuous potentials. A byproduct of our procedure is an associated lower bound 
on the maximal kissing number, which is superior to the currently best known result. By 
no means is the choice (j81J) optimal. For example, one may be able to improve our putative 
lower bound by allowing the test function to be some positive function smaller than unity for 
1 < r < a. Of course, it would be desirable to choose test functions that make asymptotic 
analysis exactly rather than numerically tractable. 

Our putative exponential improvement over Minkowksi's bound is striking and should 
provide some impetus to determine the validity of Conjecture 15.11 As a first step in this 
direction, it would be fruitful if one could show that for sufficiently small densities, the two 
standard nonnegativity conditions on the pair correlation function g 2 are sufficient to assure 
the existence of a point process, whether it is a sphere packing or not. Another problem 
worth pursuing is the demonstration of the existence of a pair interaction potential in M. d 
corresponding to a sphere packing for a given p and g 2 provided that p and g 2 are sufficiently 
small. Such a proof may be possible following the methods that Koralov |Kor05j used for the 
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lattice setting. It would also be profitable to pursue the construction of disordered sphere 
packings with densities that exceed l/2 d for sufficiently large d. 
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APPENDIX 



An asymptotic expression for J u (x) when v is large and x > v is given by |Wat58j 



cos[u u (x) - 7t/4] + O 



3x 2 + 1v 2 
12(x 2 - v 2 ) 



where 



and 



A v {x) 



Tl\/x 2 — V 2 



1/2 



(A-1) 



(A-2) 



ijj v {x) = V x 2 — v 2 — z/cos 1 {v/x). (A-3) 

The function A„(x) cos[uj u (x) — 7r/4] in ()A-ljl actually represents the dominant term in the 
asymptotic expansion (4) on page 244 of Watson |Wat58j for J v {x) when v is large and 
x > v and A v {x) multiplied by the hig-O term represents the largest absolute error when 
this dominant term is used to estimate J v (x). A problem of central concern is an estimate of 
J v {x) in the vicinity of its first positive zero xq when v is large. The first positive zero has 
the asymptotic expansion |Olv60j 

l0 = „ + ai ,V3 + I |_ + £i + o( J_ )! (A . 4) 

where 

01 = 1.8557571..., a 2 = 1.033150..., a 3 = -0.003971 ... . (A-5) 

Expanding J v (x) in a Taylor series about x = Xq and neglecting quadratic and higher-order 
terms gives the linear estimate 

JA X ) - ^[ J f-l( x o) - J v +l{xo)](x - Xq), (A-6) 
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where we take the Bessel functions on the right side to be given by the asymptotic forms 

3x 2 + 2v 2 



and 



Jv+i(x ) = A u+ i(x ) 



J v -i(x ) = A v „i(xq) 



cos[uj v+ i(xq) -7r/4} + 



COS^-i^o) - 7T/4] + O 



12(x 2 - v 2 ) 

3x 2 + 2u 2 
I2(x 2 - v 2 ) 



We will also need to consider the related functions 

1, 



J„+i(x) ~ - [Ju{yo) - Jv+2{yo)]{x - y ), 



(A-7) 



(A-8) 



(A-9) 



Jv-l{x) ^ -[J U -2(Z ) - J v (zq)](x - Zq), 



(A-10) 

where y and z are the first positive zeros of J v+i (x) and J u _i(x), respectively, which are 
asymptotically given by 



z/ + aiz/ 1/3 - 1 + 



Z/V3 3z/ 2/3 

a 2 



V 



«2 + 
V 



(A-ll) 
(A-12) 



I/V3 3Z/ 2 / 3 1/ V 4 / 37 ' 

Note that the asymptotic expressions for the zeros given here for v = 100 (d = 200) predict 
x = 108.8362071, y = 109.8641774 and z = 107.8082369, which are in excellent agreement 
with exact results x = 108.8361659, y = 109.8640469 and z = 107.8081033 as obtained 
from Maple. 

Using Maple and the results above, we obtain the following asymptotic expansions: 



1 



-[Ju-i(x ) - J v+1 (x )] 



■^.Mvo) - Ju+2{yo)\ 

\[J V -2{Z ) - J V {Z )\ 



^2/3 ^ ^4/3 ^ ' 1 > I ' 



I/ 2 / 3 1/V3 3^5/3 



Ci C2 2Ci / 1 

i/^ 



I/ 2 / 3 z/4/3 3z/ 5/3 



(A-13) 

(A-14) 
(A-15) 



where 
Ci 



2 l/4 



A/2/i(ai) + 8af7 2 (ai) 



C = 



Sy^a^ 4 

2 3 / 4 [1152a? - 3840a 4 a 2 - 180a 3 + 600aia 2 - 225] j x {a x ) + 2 1 / 4 [3072ai /2 - 200a 3 / 2 



(A-16) 



. (2 ai ) 3 /2\ /(2ai 
/i(ai) = sin ( J ) + cos 



3/2 



3840v^Faf /4 



/ 2 (ai) = sin 



(2a x 



1 3/2 



COS 



(2a x 



,3/2 



(A-17) 
(A-18) 
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Thus, substitution of the values for the constants a\ and a 2 into the expressions above yields 
the estimates 

d = -1.104938082, C 2 = 1.627074727. (A-19) 

For d = 200 [y = 100), for example, the asymptotic expansions (A-13) - (A-15) predict 
-0.04778125640, -0.04743934518 and -0.04812316762, respectively. This is to compared 
to the corresponding exact results: -0.04829366129, -0.04799533693 and -0.04859672879. 
Note that although the estimates of the constants C\ and C 2 given by (|A-19J) involve a small 
error due to our use of only the dominant asymptotic term (jA-l|) . the functional forms of the 
asymptotic expansions (A-13) - (A-15) are exact. 

We can show that the exact expressions for the constants C\ and C 2 appearing in (A-13) 
- (A-15) are rapidly converging asymptotic expansions in inverse powers of the constant a\ 
appearing in (|A-4jl . For example, using expansion (4) on page 244 of Watson, we find that 
the first three terms of the expansion of C\ is given by 



Ci 



C 11 C*12 C 1 3 m [ 1 



where 



2 l/4 



5/4 



x/2/ 1 (a 1 ) + 8^ /2 / 2 (ai 



11/4 



17/4 



23/4 



(A-20) 



5 • 2 1 / 4 



C 



4 x /2/ 1 (a 1 )-7a? /2 / 2 (ai 



7T 



12 



384^ 



385 ■ 2 1 / 4 



C 



13>/2/i(ai) + 8al /2 f 2 ( ai ) 



13 



221184v^F 

The first term of the expansion (jA-20|) is the dominant one and is identical to the estimate 
given in (jA-16|) . The first term of (jA-20|) is about 66 times larger than the second term in 
absolute value and the second term is about 7 times larger than the third term in absolute 
value. Substitution of the constant a\ into (jA-20|) yields the more refined estimate 



Ci 



-1.123958144. 



(A-21) 



This refined estimate differs from the dominant first-term estimate given in ()A-19j) in the third 
significant figure. One could continue correcting this estimate by including additional terms 
in the asymptotic expansion but this quickly becomes tedious and is not necessary because, 
as we show at the end of Section 5, the precise value of C\ is not relevant for the putative 
exponential improvement of Minkowski's lower bound on the density, as specified by (|116|) . 
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